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ABSTRACT 

A previous paper [1] presented an approximate 
analytical treatment of this problem. That treat- 
ment resulted in an expression for the apparent 
absorptance of the condensate that compared well 
with experimental data reported elsewhere. The ex- 
pression snowed Ihat under certain enecified condi- 
tions, the presence of a condensate could alter the 
heat transported to a surface significantly, when 
compared with the heat transported in the absence of 
a condensate. In obtaining a closed form solution, 
however, assumptions were made that might have had 
the effect of limiting the solution. 

The present paper refines and extend© the 
earlier analysis by: 

(1) allowing for an explicitly varying conden- 

sate depth 

(2) considering a range of substrate tempera- 

tures 

(3) treating emission frem the condensate more 

precisely 

In addition, a variety of boundary conditions have 
been imposed at the condensate -substrate interface 
in order to extend the applicability of the anal- 
ysis. 


The results to date indicate that the previous 
assumptions did not invalidate the earlier analysis 
for the broader class of problems considered here. 
Even at a substrate temperature of 250° K, the rate 
of sublimation from the condensate was so small that 
the variable condensate depth had a negligible 
effect upon the temperature distribution, as also 
did the emission. The effect of higher substrate 
temperature was to increase the heat loss coefficient 
h of the condensate. Problems with other boundary 
conditions are being investigated. 


INTRODUCTION 

A previous paper [1] presented an approximate 
analytical treatment of the transport of thermal 
energy in an H 2 0 deposit exposed to thermal radia- 
tion in a vacuum. The deposit (frost) was assumed 
to have grown on the exterior surface of a cryogenic 
fuel tank. The analysis showed that the presence of 
an H 2 0 deposit could increase the heat transported 
through a reflective surface by as much as a factor 
of eight over the value transported through the hare 
surface. To achieve this closed form solution, 
however, two principle assumptions were made: (l) 
the change in deposit thickness due to sublimation 
was assumed to be negligible, so that the differ- 
ential equation of heating showed no explicit de- 
pendence upon the rate of sublimation, and (2) the 
radiant energy emitted by the frost was combined 
with the heat lost by sublimation and both were 
treated in the boundary conditions. The purpose of 
the present paper is to refine the analysis by re- 
moving these assumptions, and further, to apply the 
refined analysis to the case of cryode posits formed 
upon shadow shields and within multifoil insulation 
systems . 


ANALYSIS 

The situation to be analyzed is presented in 
figure 1. Thermal »■*.- liation is incident upon a 
nonopaque solid, the cryodeposit, which covers ihe 
exterior metallic surface of a cryogenic reservoir. 
Of the radiant energy incident upon the frost, cer- 
tain portions will he 

(1) Reflected to space at the vacuum-frost 

interface 

(2) Consumed in sublimation of the frost 

(3) Absorbed in the frost raising the frost 

temperature 

(4) Absorbed in the frost and conducted to the 

reservoir 

(5) Absorbed in the frost and reradiated to 

space 

(6) Transmitted through the frost directly to 

the metallic substrate 

Of the radiation transmitted through the frost di- 
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rectly to the substrate, a portion will be reflected 
back to the frost again as determined by the reflec- 
tance at the frost -metal interface. By considering 
the energy balance for a differential volume element 
(shown in fig. 1) the following heat transfer equa- 
tion is obtained: 



' a 2 iA 

.ay 2 / 


^ - Pl2) w ine K a 


X 


-K [H-s-y] 


-K [H-s+y] 


+ p. 


23 


- K e o(T 4 - T*j (1) 

where K a is the absorption coefficient averaged 
over a spectral region of interest (the "box" 
approximation of [2]) and W{ nc is the incident 
radiant energy in that spectral region. This defi- 
nition of K a and W.[ nc is well suited to ma- 
terials that are absorbent over limited spectral 
regions. In such cases, the contribution of each 
absorption band would simply be added. The absorp- 
tion coefficient, #c e , is also averaged, but over the 
spectral region of emission. In this analysis, /< a , 
K e , and all other physico-chemical properties of the 
deposit are considered to be known. Symbols are de- 
fined in the appendix and physical properties used 
in the calculations are presented in table I. 


TABLE I. - PHYSICAL PROPERTIES OF WATER CRYODE POSIT 


Property 

Value 

Density of ice 1 
Heat of sublimation 1 
Specific heat 2 
Thermal conductivity 3 
Absorption coefficient for 
solar radiation 4 
Absorption coefficient for 
290° K blackbody 4 
Absorption coefficient for 
77° K blackbody 4 
Emission for 250° K case 

0.917 g/cu cm 
3050 j/g 
0.697 j/(g)(°K) 

0.078 j/(cm)(sec)(°K) 
0. 0305 /q 

0. 1578 /fi 

0.0640 /u 

0. 1578/|i 


W. 3. 

Ref. 4 
^Ref . 5 
4 Ref. 6 


In physical terms, the boundary conditions on 
the solution of the problem are: 


The analytical statement of these conditions is: 


at y = 0 

T = T c 

at y = H - s 

•Kg) ■ 

at t = 0 

II 

J-3 

O 


Equation (l) is nondimens ionali zed by means of the 
following substitutions 

T - T p kkot 

q = . | = K a [ (H - s) - y]; t = ; 

i c M p 

ff s (t) s * a s(t); h s Kq H 


Thus, 0 is the temperature rise scaled by the sub- 
strate temperature, £ is the optical depth of the 
frost, and t is the Fourier number of classical 
conduction heat transfer. The quantities a g and 
h are introduced for convenience. In terms of 
these new variables, the heat transfer equation is 

de_. /d0\ /d 2 e\ f . -[2(h-a s )-|]l 

57 =<,s W/ = (^j + Q | e ' l+p 23 e 8 ]-3‘ e 

(la) 

where 

(1 - p 12 )W inc _ 

' 9 “kK 2 T 

a c 

In order to show the extensions involved in the 
present analysis, the present heat transfer equa- 
tion (la) is compared with the equivalent equation 
solved in [1]. In terms of the present symbols 
that equation is 

d0 /% 2 o\ r -(2h-|)l 

37 \d?j + T + P 23 6 J 

This last equation is obtainable from (la) if: 

(1) The contribution of the emission term Q* 
is not included in the equation. 

(2) a s « h 

(3) & s(|f) «(|f) 


(1) The substrate temperature is constant and 
equal to T Q 

(2) The frost may lose heat and mass at the 
vacuum surface as a result of sublimation. 

(3) The frost is in a quiescent steady-state 
condition at a definite thickness H and definite 
temperature T c prior to the incidence of thermal 
and radiation. 


These conditions were assumed in [l]; however, these 
assumptions are not made in the present analysis. 

The partial differential equation that is 
solved in the present analysis is 

- & S (H) = + { e_l + P 23 e- [2(h -° s) - n |- 

(2) 


\ 
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where \|r = (e/Q) is introduced for convenience. As 
shown in detail in [l], the houndary conditions may 
be linearized to result in the following: 



at |=h-cr s i = 0 

at t = 0 1 ^ = 0 

Finally , there is a subsidiary relation for a g 



Several analytical solutions of (2) were in- 
vestigated in a cursory fashion. Notable among 
these was a solution in terms of a transformation 
first suggested by Paterson [7] and more recently by 
Sneddon [8]. By means of the transformation 

u ^ \jre" H ^' t 

equation (2) becomes 



S 2 u 

se 2 



+ p 23 6 


-[2(h-a )-£] 


As Q 1 » 0[10" 5 ], the solution to this equation 
would be expected to be very similar to that pre- 
sented in [1] were it not for the convective term, 

, whose contribution could be significant. 


No exact solution was found fur the equation contain- 
ing the convective term; therefore , a numerical so- 
lution w«s sought. In this way. a measure of the 
contribution of each term could be obtained. 


Figure 2 presents the spatial distribution of 
the dimensionless temperature rise through a water 
frost on a 250° K substrate. The physical proper- 
ties of the frost are listed in table I. The reason 
for this choice of substrate temperature will be 
clarified subsequently. The Hue at which ste^y- 
state was achieved was arbitrarily defined as that 
time after which there was no change in the third 
significant figure of v|f. In figure 2, the abscissa 
is £/(h - a s ); however, for all cases considered, 
it was found that a g < 10 ~ 2 h. 

The complete energy flux to the reservoir is 
given by 

CL = (1 - PoJ(W. - W* ) + cl Jj] (1 - p, J 
tot ' 23' v inc me tot me' 12 


In a particular case, the dominant absorption term 
will depend upon the spectral properties of the 
deposit ( i. e . , k as a function of wavelength). The 
expression for t is 


tot 


-it), 


^=h-a s 


(1 


K 23 


)e 


-Ch-a ) 


Thus, the curves for a^.^, presented in figure 3, 
give the heat flux to the substrate by conduction 
fron the frost and radiant transmission through the 
frost in the spectral region appropriate for w! nc * 
These curves, which were first presented in [9], 
give the heat transfer at steady state as defined 
previously. In the present paper, a rather high 
value of substrate temperature (250° K) was chosen 
in order to assess the effect of an appreciable 
mass-transfer rate. Even for this case, a s « h. 

As an example for h = 0.125 with the deposit at 
250° K subjected to 290° K blackbody radiation, 
cf s s 0.15x10' 7 at steady- state. Thus, the assump- 
tion cr s « h is valid for this case. Recall that 
[1], based upon that assumption, did achieve a 
truly steady-state solution. Therefore, the results 
presented here tend to confirm the results pre- 
sented in [l], and further, the steady-state results 
as here defined are equivalent to the steady- state 
results of [l]. For the H 2 0 deposit considered 
here, the wavelength region of interest was taken 
to be 2 to 20 microns for absorption and 2 to 100 
microns for emission from a substrate whose temper- 
ature ranged from 20° to 250° K. The spectral 
regions were chosen to include the principle ab- 
sorption bands of the deposit. A temperature of 
250° K could be obtained by the bulkhead of a 
spacecraft on the night side of an Earth orbit or 
by an interior surface of a shadow shield array. 
Thus, these findings indicate that a deposit, once 
formed, could remain stable on such surfaces. 

The results of a more complete analysis tend 
to confirm the validity of the assumptions made in 
an earlier analysis, namely, that o g « h, and 

that the effects of Q' and c s are negligible. 

Several applications of cryodeposits now will be 
investigated. 


APPLICATIONS 

In the case of bare surface reflectance shield- 
ing, the effect of h »-vy.-.depocit is shown in fig- 
ure 3. Fran the analysis of [1] and figure 3, the 
total heat transport due to incident energy in the 
range W± nc decreases as l/h for h > 1.0 and 
h » 1.0. On this basis, [9] suggested that a 
deposit could be used to reduce heat transport to a 
surface that had been degraded. To accomplish 
+>n's. the minimum deposit optical thickness is 
given by the condition that 

(1 - P 12 ) a tot < 1 ’ p 23 

Bare surfaces of cryore servo ir s , however, are 
rarely, if ever, exposed to direct solar radiation. 
Most often, at least one additional surface exists. 
Depending upon the temperature of that surface, the 
radiation it emits could be in just that spectral 
region for which the frost is most absorbent (e.g., 
2 to 20 |i in the case of H 2 0). Specific examples 
have been calculated for the configuration of 
shields and cryodeposits shown in figure 4. The 
steady- state thermal energy flux through a two 
shield array is obtained from the simultaneous 
solution cf the following equations: 

aW inc ‘ *1^1 = E 12 F 12 a ( T l ‘ T 2^ + £ jS X " *12^1 


3 



E 12 F 12 a ( T l " T 2^ ' e 2^ ' p 12) oT 2 

" E 2E F 2H 0(T t - T R> + - F 2rK 


where 



and where a and W inc refer to the complete solar 
spectrum. The present calculations are based upon 
Tr = 20° K, F]_2 = 0.382, and € = 0.045 for the hare 
shield surface (e.g. , aluminum foil). It is further 
assumed that the surfaces with deposits have emit- 
tances of 0.090. As shown in figure 3, this emit- 
tance would correspond to a frost thickness of less 
than 3(j for K a = 0.1578 per p. The resultant 
heat fluxes to the cryoreservoir are shown in fig- 
ure 4 for each configuration. In both cryodeposit 
cases, the heat flux is increased by over 30 per- 
cent due to the presence of the deposit. It is 
acknowledged that the figures quoted are based upon 
a steady-state calculation. Neverthless, the tem- 
perature attained by the inner shield and the tem- 
perature of the reservoir are such that a deposit, 
once formed, can remain stable. 

The effect of the boundary condition of 
constant -reservoir temperature was investigated to 
determine the applicability of the present analysis 
(i.e., the solution of eq. (2)) to the shielding 
problems. As an extreme, the rate of sublimation 
was calculated for the case in which the inner frost 
surface was insulated against conduction heat trans- 
fer. Even in this case, a steady-state temperature 
distribution was obtained while a & « h. Thus, the 
calculations are valid for the case of shadow 
shields. 


not appear that the condensate would adversely in- 
fluence the effectiveness of a multilayer insulation 
system. 

In conclusion, a numerical solution of the 
complete heat balance equation has substantiated an 
earlier preliminary analysis of heat transport in a 
nonopaque frost layer. In particular, sublimation 
was found to be sufficiently small for all cases 
considered such that o s « h. Thus, it was con- 
cluded that the previous analysis was adequate. 

With this in mind, the results of [9] were reviewed 
and extended to the case of multilayer insulators. 

It was found that the effectiveness of the multi- 
layer insulation system probably would not be ad- 
versely affected by the presence of a cryogenic 
condensate. 


APPENDIX - SYMBOLS 

0 ^ specific heat of deposit 

radiant exchange factor between parallel 
plane surface i, j 

configuration factor for parallel plane 
J circular surfaces ij 

H initial frost thickness 

h initial frost optical thickness /c & H 

h linearized heat loss coefficient 

k thermal conductivity of deposit 

L g heat of sublimation of deposit 

m rate of deposit mass loss due to sublima- 

tion 


P. J. Perkins, of NASA Lewis, suggests that 
multilayer insulation be vacuum sealed and filled 
with a readily condensable gas such as COg in order 
to provide a cryopumped vacuum within the insulation 
while in the Earth's atmosphere. In order to inves- 
tigate the effect of the condensate upon the effec- 
tiveness of the system, a simplified multilayer in- 
sulator was investigated. It was assumed that five 
shielding layers separated by vacuum spacing were 
used to shield a fuel tank. In the steady-state 
the heat flux to the reservoir is given by 


Q dimensionless incident radiant energy 

(^. ot complete energy flux to reservoir 

Q* dimensionless emitted radiant energy 

s sublimed deposit depth 

s rate of sublimation, real time 

T temperature 


£lOT* 


T c temperature of substrate 


Again, for bare surfaces € = 0.045; for frost 
covered surfaces, e = 0.090. Since W net is also 
tbe net flux between neighboring foils, there re- 
sults five equations for five unknown foil tempera- 
tures (T r is presumed known). For the case of 
bare shields, w net ~ 99x10" 3 watts per square cen- 

timeter and for the case where a condensate forms on 
the reservoir and the innermost shield, 

W ne £ » 2. 18x10" 3 watts per square centimeter, a 
10 -percent increase in the hoiloff heat flux. The 
model considered here does not allow for spacer con- 
ductance or gas conductance between shields whose 
temperatures would preclude condensation. Both of 
these conductances could easily increase the heat 
flux by more than 10 percent. Consequently, it does 


Tr, reservoir temperature 

.K 

t time 

w inc total incident radiant energy 

W inc radiant energy in spectral region for which 

frost absorbs 

y spatial variable 

a absorptance 

a. . total dimensionless heat flux 

tot 

e emittance 


1 


4 


E-34S3 



E-342 




A 

S 


P 



T 


* 


dimensionless temperature rise 

absorption coefficient, averaged 
spectral region of interest 



over 


See [2], 


the 1v box” approximation. 

absorption coefficient, average over spec- 
tral region appropriate to the emitted 
energy 

wavelength 

dimensionless spatial variable 

density of deposit 

reflectance at exterior interface 

reflectance at interior interface 

Stefan-Boltzmann constant 

optical depth loss due to sublimation 

rate of sublimation, transformed time 

transformed time (Fourier number) 

dependent variable for dimensionless tem- 
perature rise 
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Figure 1. - Cryodeposit and substrate complex. 



Fig. 2 . - Steady- state distribution of dimensionless tem- 
perature rise. 
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Fig. 3. - Total dimensionless heat transport through a 250° K substrate for solar and 290° K 
blackbody radiation. 
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Case 1 . - No cryodeposits present. 
Incident flux, 0. 14W/cm 2 ; solar absorp- 
tance, 0. 15 for first shield; emittance, 
0.045; net boiloff flux, 2. 68x10-4 W/cm 2 ; 
middle shield temperature, T 2 , 272° K. 


Case 2 . - Deposit on inner shield 
and reservoir. Assumed emittance of 
deposit coated surfaces, 0.090; net boil- 
off flux, 3. 52x10' 4 W /cm 2 ; middle shield 
temperature, T 2 , 246° K. 


Cryodeposit 


N \ 

' | Cryogenic 
? reservoir 


Case 3 . - Deposit only on reservoir. 
Assumed reservoir emittance, 0.090; 
net boiloff flux, 3. 50xl0 -4 W/cm 2 ; middle 
shield temperature, T 2 , 268° K. 

Fig. 4. - Effect of cryodeposition on shadow shields. 
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